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Abstract 


Protection of astronauts and instrumentation from galactic cosmic rays (GCR) 
and solar particle events (SPE) in the harsh environment of space is of prime 
importance in the design of personal shielding, spacecraft, and mission plan- 
ning. Early entry of radiation constraints into the design process enables op- 
timal shielding strategies, hut demands efficient and accurate tools that can 
he used by design engineers in every phase of an evolving space project. The 
radiation transport code, HZETRN, is an efficient tool for analyzing the shield- 
ing effectiveness of materials exposed to space radiation. A systematic effort 
of verification and validation is now underway to guantify the accuracy of 
HZETRN and elucidate any weaknesses. In this paper, HZETRN is compared 
to the Monte Carlo codes HETC-HEDS and ELUKA, for a shield/target con- 
figuration comprised of a 20 gfcmf Aluminum slab in front of a 30 gfcmf 
slab of water exposed to the Eebruary 1956 SPE, as modeled by the Web- 
ber spectrum. Neutron and proton fluence spectra, as well as dose and dose 
eguivalent values, are compared at various depths in the water target. Elu- 
ence spectra produced by HZETRN are also given for ^H, ^H, ^He and ^He 
ions at various depths in the water target. This study shows that there are 
many regions where HZETRN agrees with both HETC-HEDS and ELUKA for 
this shield/target configuration and the SPE environment. However, there are 
also regions where there are appreciable differences between the three computer 
codes. 


1 Introduction 

One of the objectives of the NASA Strategic Plan for the human exploration of space 
is to develop space radiation protection for both astronauts and instrumentation [1], An 
important step in achieving this goal is to provide design engineers with an efficient 
radiation transport code with an estimate of accuracy under a variety of design scenarios. 
Such a code can then be used in every phase of vehicle, shelter, and mission design to 
enable optimal shielding strategies that satisfy radiation constraints and objectives. 

The HZETRN (High charge(Z) and Energy TRaNsport) computer code has been used 
for radiation analysis under a variety of shielding conditions in solar particle event (SPE), 
galactic cosmic ray (GCR), and low earth orbit (LEO) environments. While the code has 
endured several rounds of verihcation and validation in these environments [2, 3, 4], most 
of the comparisons focused on integrated quantities, such as dose or dose equivalent, and 
individual ion fluences were not examined in detail. Though dose or dose equivalent were 
generally viewed as sufficient tests for evaluating code accuracy, recent interest in fluence 
based approaches to radiation risk assessment demand a higher degree of accuracy from 
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the code [5]. This, along with recent improvements to some of the underlying transport 
models and numerical procedures [6, 7, 8], provides an opportunity to expose HZETRN 
to another round of verihcation and validation benchmarks. 

In this paper, a 20 g/cm^ Aluminum slab shield and a 30 g/cm^ water slab target 
exposed to the February 1956 SPE spectrum, as modeled by Webber [9, 10] is used as a 
benchmark case to compare HZETRN to the Monte Carlo codes HETC - HEDS (High 
Energy Transport Code - Human Exploration and Development of Space) [11, 12] and 
FLUKA [13, 14]. The 20 g/cm^ Aluminum slab is a common thickness for space vehicles 
and the 30 g/cm^ water slab simulates body tissue. 

While HZETRN has historically transported all particles in the straight ahead direc- 
tion, recent efforts have provided a bi-directional neutron transport model that is fully 
coupled to the existing transport algorithm. References [6, 8, 15, 16] provide a complete 
description of the straight ahead transport procedure, and a brief description of the the 
bi-directional neutron transport model and coupling mechanisms used by HZETRN is 
provided. An overview of the Monte Carlo codes is also given with suitable references to 
provide more detailed information. Forward and backward neutron fluences and proton 
fluences are generated using all three codes and results are compared at various depths in 
the water target. Dose and dose equivalent comparisons between HZETRN, HETC-HEDS 
and FLUKA are given. Since HZETRN also transports ^H, ^H, ^He and ^He for any SPE 
environment, these fluence spectra are also presented. 

2 HZETRN 

HZETRN is a suite of codes containing a numerical solution for the Boltzmann trans- 
port equation with the continuous slowing down and straight ahead approximations. The 
code has evolved over a 25 year period from the work of John W. Wilson and co-workers at 
the NASA Langley Research Center [15]. Several corrections and modihcations have been 
added in recent years. A more robust method for handling neutron elastic interactions 
has been included [6], and a bi-directional neutron transport model has been fully cou- 
pled into the code. Light ion (A < 5) cross sections have been extracted to allow greater 
computational flexibility and efficiency. These cross sections were previously imbedded in 
the code. 

The Boltzmann transport equation, with the continuous slowing down and straight 
ahead approximations, is given by [15] 

B[4>^{x,E)] = J2 [ a,k{E,E')M^,E')dE' , ( 1 ) 

k 


2 



where 


B[4>j{x,E)] 


A 

dx 




S,{E) + a,{E) 


(t)j{x,E) 


( 2 ) 


is a differential operator. Here, (f)j{x, E) denotes the fluence of type j particles at depth x 
having kinetic energy E, Aj is the atomic mass of a type j particle, Sj{E) is the stopping 
power of the jth particle with kinetic energy E, <Jj{E) is the macroscopic cross section, 
and ajk{E,E') are the production cross sections for interactions where a type k particle 
with energy E' produces a type j particle with kinetic energy E. 

The following are some notes concerning the calculation of the numerical solutions to 
equation (1). For heavy ions {A > 4), it is noted that projectile fragments, after leaving 
a collision event, will have energy very near that of the projectile. Therefore, spectral 
distributions in the production cross sections become very narrow and allow for easy 
evaluation of the source integral on the right hand side of equation (1) [15]. Further, target 
fragments created by heavy ion projectiles are not explicitly transported due to their low 
range (they are accounted for in dose calculations) and so the summation in equation (1) 
is taken over k such that Ak > Aj [15, 17]. For light ions, neither of these approximations 
is valid. The broad energy distribution in collision events prohibits simplihcation of the 
source integral, and both projectile and target fragments are transported, so that the 
summation over projectile particle type must be taken over all light ions [17]. Since this 
paper deals with the February 1956 Webber SPE spectrum, for which there is a negligible 
heavy ion component, no more will be said about heavy ion transport in HZETRN. 

Wilson et ah have developed an accurate marching algorithm to transport light ions 
in the straight ahead direction [8, 15, 16]. The formulation and implementation of the 
procedure has not changed much in recent years, but Slaba et al. [6] have provided a 
more robust method of handling low energy neutron elastic interactions, and light ion 
cross sections have been successfully extracted from the code to allow greater computa- 
tional efficiency and most importantly, the ability to interchange cross section models and 
databases as necessary to carry out future verihcation and validation efforts. 

The most notable modihcation to the HZETRN code came in recent years with the 
culmination of several works [18] - [24]. Slaba et al. [7] have developed an efficient 
and accurate bi-directional neutron transport model that is completely coupled into the 
existing light ion transport procedure. The fluence (f)j is split into a straight ahead (sa) 
and isotropic (iso) component by writing 


0,(o;,E) = 0“°(a;,E) + 0f(a;,E), (3) 

where it has been demonstrated that the straight ahead component is associated with 
higher energy primary or cascade projectile fragments produced in nearly the same di- 
rection of the projectile, and the isotropic component is associated with particles ejected 
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from excited projectile fragments or evaporation. Similarly, the neutron nuclear reactive 
production cross section is split into straight ahead and isotropic components as 


ank{E,E') 


= (J. 


nk 




( 4 ) 


where the superscript r denotes the nuclear reactive component of the neutron production 
cross section and the subscript n denotes neutrons. The straight ahead component of the 
particle fluences are first computed by solving 

B[<t^f] = T.f E') dE' , (5) 

k 


using the light ion marching algorithm developed by Wilson et ah [8, 15, 16], with 
Kjk{E,E') given by the piecewise dehnition 


Kjk{E,E') 


^:ik{E,E'), j=n 

ajk{E,E'), j^n. 


( 6 ) 


The isotropic neutron source term, 

E) = J2[ oT‘(E, E'W,-(x, E') dE' , ( 7 ) 

T-'® 

is then calculated, and the isotropic neutron fluence is decomposed into a forward 
and backward component by writing 

^:°(x,E) = 4>’n(x,E)+^Hx,E) ■ ( 8 ) 

A collocation method [7, 23, 24] is hnally used to obtain a Neumann series solution to the 
directionally coupled neutron transport model, 

= / a^^^n{E,E')<t>i{x,E')dE' 

J E 

+ [ E')(j)’’^{x, E') dE' + E) , (9) 

J E 

= / ‘x<;»(E,E')4,i{x,E’)dE' 

JE 

+ l^alc„\E,E')<pHx,E')dE' + . ( 10 ) 

The production cross section, a^\E, E'), represents nuclear reactive and elastic interac- 
tions, in which the scattering angle between the pre - collision and post-collision neutron 
is between [0,7t/2] with respect to the propagation direction of the pre-collision neutron. 


dx 


+ O’n(E) 


(l^i{x,E) 


d . 


d > i ( x , E ) 
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Conversely, a^\E,E') represents nuclear reactive and elastic interactions in which the 
scattering angle between the pre-collision and post-collision neutron is between [7r/2,7r] 
with respect to the propagation direction of the pre-collision neutron. This implies that 
the second integral on the right hand side of equation (9) represents forward propagating 
neutrons produced by backward propagating neutrons, and the second integral on the 
right hand side of equation (10) represents backward propagating neutrons produced by 
forward propagating neutrons. In order to couple this bi-directional neutron solution back 
into the light ion transport procedure, the isotropic light ion source terms. 


vrix, E)= [ a,^{E, E') dE' , (11) 

J E 

are calculated, and the isotropic component of the light ion fluences are obtained by 
computing 


<i>r{x.E) = 


P,{E,)S,{E,) 
P,{E)S,{E) ^ 

E. A,P,{E') 


+ 


Ie Pj{E)Sj{E) 
where Rj{E) is given by the range-energy relation. 


7]j{x + Rj{E) - Rj{E'), E') dE' , 


( 12 ) 


R,{E) = Ai / 


E dE' 


Sj{E') ’ 

and Pj{E) is the nuclear survival probability dehned by 


(13) 


Pj{E) = exp 


-A- r^mdE' 


lo S,{E') 


(14) 


with the residual energy given by E^ = R~^[x + Rj{E)]. More detailed information 
regarding these transport models and numerical implementations can be found in the 
references [7, 18]. 


3 Monte Carlo Codes 

The term, Monte Carlo, is used to describe a random walk algorithm that simulates an 
event and then performs a statistical analysis of the results. Monte Carlo techniques are 
typically constructed whenever there is a high dimensional integral equation that cannot 
be solved either analytically or numerically with deterministic methods. To reduce the 
computing time involved in Monte Carlo techniques, variance reduction is often employed. 
A brief description of Monte Carlo methods applied to radiation transport can be found 
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in reference [25] . In this paper, the results from the two Monte Carlo codes HETC-HEDS 
[11, 12] and FLUKA [13, 14], are used to verify the results calculated by the deterministic 
code HZETRN. The following is a brief description of HETC-HEDS and FLUKA used in 
this study. 

3.1 HETC-HEDS 

The HETC-HEDS computer code is a Monte Carlo based solution method designed 
specihcally for solving space radiation problems, [11, 12] associated with secondary particle 
helds produced by space radiation interacting with shielding and equipment. It is a 
three dimensional generalized radiation transport code capable of handling and analyzing 
radiation helds which affect critical body organs of astronauts such as bone marrow and 
the central nervous system. 

HETC-HEDS can be applied to a wide range of particle species and energies. It con- 
tains a heavy ion collision event generator capable of tracking nuclear interactions and 
performing statistical analysis of the data. It simulates particle interactions by using a 
pseudo random number generator, along with the appropriate physics, to follow trajec- 
tories of primary particles and all secondary particles involved in the nuclear collision of 
galactic cosmic rays and solar event particles interacting with shielding material, biological 
organisms, and electronic equipment. 

The geometry input is simple combinatorial volumes. The cross sections are gener- 
ated internally, as part of the software. The HETC-HEDS computer code employs all 
particles of interest for space radiation. In particular, HETC-HEDS considers interac- 
tions of protons, neutrons, n~ , fx~ , light ions and heavy ions. These particles 

can be arbitrarily assigned position, angle, and energy throughout a spatial boundary. 
In essence, this Monte Carlo code follows each particle in a cascade until it undergoes a 
nuclear collision, absorption, decays, escapes from the spatial boundary, or is eliminated 
by crossing a domain variable cutoff. The nuclear reactions and processes are accounted 
for by using appropriate physical models to handle such things as energy losses, range 
straggling. Coulomb scattering, etc. Both elastic and nonelastic collisions are considered 
using energy and nucleon conservation principles. 

This Monte Carlo package has no capability to determine when the solution method 
has converged on the answer sought within some conhdence interval. It does provide the 
means for predicting the interaction product yields, production angles, and energies using 
nuclear models for transport processes. It has been extensively used for code verihcation, 
bench marking and testing against available laboratory beam data. 

3.2 FLUKA 

The FLUKA computer code is a general purpose Monte Carlo computer program used 
for calculating particle transport and interaction with various materials. It has the ability 
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to transport and interact all elementary hadrons, light and heavy ions, and electrons and 
photons over an energy range which extends up to 10^ TeV for all particles, and down 
to thermal energies for neutrons [13, 14], The code has built in capabilities for scoring 
particle fluences, yields, and energy deposition over arbitrary three dimensional meshes, 
both on an event by event basis or averaged over a large number of histories. 

This code has been extensively benchmarked against available accelerator and cosmic 
ray experimental data, at beam energies as low as a few MeV and as large as cosmic ray 
energies. These spectra can be modulated within FLUKA according to an arbitrary solar 
activity modulation parameter, or, for dates in the past, using the actual solar activity 
as measured by ground based neutron counters. 

Three different Earth magnetic field descriptions of increasing complexity can be used. 
Geomagnetic cut offs can be input or calculated. Spectra representative of some of the 
largest solar particle events are also pre-built into the code. In addition to galactic cosmic 
rays and solar particle events, this code has a wide range of particle physics applica- 
tion areas where it can be employed. It has demonstrated that it can simulate nuclear 
interactions and events with great accuracy in several helds, including particle physics, 
dosimetry, accelerator driven nuclear systems, detector characterization, aircraft crew 
dosimetry, shielding, and hadron therapy. The code can simulate interactions involving 
over sixty different particle types and can handle complex geometries. The geometry in- 
put has two levels: the basic one consisting of combinatorial volumes, and a second level, 
the so called “lattice” level, which allows repetition and arbitrary spatial placement of 
the complex object dehned at the first level. It also allows for a voxel (a volume element, 
representing a value on a regular grid in three dimensional space) input stream which 
can be combined with a standard combinatorial geometry input, allowing for example to 
embed a detailed voxel description of a human being, such as derived from a computed 
tomography (CT) scan, within an arbitrarily complicated spacecraft shape. 

The nuclear models and associated cross sections are hard coded into the software 
for all particles and energies, with the exception of neutrons below 20 MeV, where a 260 
neutron group library obtained from standard international evaluated data hies is used. 
The FLUKA computer code contains many variance reduction algorithms and covers all 
particles of interest for space radiation research. For most applications, the FLUKA 
package requires no additional programming. 

The code provides a large variety of statistical techniques for analyzing nuclear inter- 
actions. The offline analysis tools provided with the code enable one to perform simple 
statistical analysis of the results. The code comes with a powerful Graphical User Inter- 
face, FLAIR, which allows an easy setup of the input stream and the geometry. It also 
provides user friendly tools for running the code and analyzing and plotting the results. 
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4 Results and Discussions 


The 1956 Webber differential and integral SPE spectrum, with 100 MV rigidity, is 
illustrated in ffgure 1. The Webber spectrum is given by 


d<p 

dE ^ 

4 > = 


10^ ^i(E+i876) + 1876)] /lOO} 

10^ exp { [239.1 - ^E{E + 1876)] /lOO} 


# part ides 
cm^ MeV ’ 

T^partides 
cm^ ’ 


( 15 ) 


and is exclusively a proton spectrum [9, 10]. 

Figures 3-19 present selected results generated by the three radiation codes used to 
transport the February, 1956 Webber SPE through a 20 g/cw? Aluminum slab shield and 
a 30 g/cw? water slab target. Figures 20 - 23 illustrate ^H, ^H, ^He and ^He ffuences 
predicted by HZETRN. These four light ions along with neutron and protons are the 
six light ions transported by HZETRN for SPE environments. They have been included 
in this paper to give a complete picture of the particles used by HZETRN to calculate 
dosimetric quantities for SPE environments. 

Observe that in many sections of the graphs presented, the data from all three com- 
puter codes are overlapping and consequently one cannot distinguish some of the markers. 
A cursory examination of the results presented in figures 3-19 allows one to conclude 
that there is a good qualitative agreement between the transport codes in many regions. 
However, there are some regions where there are significant differences between the codes, 
and the following discussion will address these issues. 

For a quantitative assessment of the dose and dose equivalent curves, we select a root 
mean square analysis of the data. Recall that the root mean square between data sets 
{Yi} and {yi] for i = 1, . . . , V is given by 


RMS 


N 


(16) 


where N is the number of ordinates in the data sets with the abscissa values being the 
same for each data set. We make a RMS comparison between the HZETRN, FLUKA and 
HETC-HEDS data sets and obtain 


RMSDose(HZETRN-HETC) 0.176 cGy , 

RMSDose(HZETRN-FLUKA) ~ 0.212 cGy , 

RMSDose(pLUKA-HETC) = cGy . (17) 


An examination of the dose equivalent curves in ffgure 4 shows there is a significant 
difference between the calculated results. We again use a root mean square analysis for a 



quantitative measure of agreement between the data sets and obtain the value 


R-MSDoseEq(HZBTRN-HETC) ~ 1-085 cSv , 

RMSDoseEq(HZBTRN-FLUKA) ~ 0.184 cSv , 

RMSDoseEq(pLUKA-HBTC) ~ 1-008 cSv . (18) 

Note that a small RMS value implies better overall agreement between the data sets. For 
the dose curves, there is good agreement between FLUKA and HETC-HEDS. For the 
dose equivalent curves, the best agreement is between HZETRN and FLUKA. Both the 
HZETRN and FLUKA computer codes use the ICRP-60 standard for calculating dose 
and dose equivalent, while HETC-HEDS uses the ICRP-25 standard. It should also be 
emphasized that HZETRN does not transport certain particles such as pions, muons, 
positrons, electrons, and photons. These particles are used in calculating dose and dose 
equivalent by HETC-HEDS and FLUKA, but not HZETRN. The contribution of these 
particles to dose and dose equivalent values can be signihcant. Also, the HETC-HEDS 
uses the NCRP 132 [26] and 142 [27] conversion factors in combination with the DABL69 
KERMA factors [28] . The HZETRN code does not use KERMA factors in the calculation 
of dose and dose equivalent. The number and type of particles used and the algorithm 
employed for calculating the dose and dose equivalent could account for the differences 
in the calculation of the dose and dose equivalent curves [25]. In calculation of the dose 
equivalent, the HZETRN code uses the lineal energy quality factor Q{L) as dehned in 
reference [29] 


Q{L) = 1, L<10keV//rm, 

Q{L) = 0.32L — 2.2, 10 keV//rm < L < 100 keV//rm , 

Q{L) = 300L“°'^, L > 100 keV//im . (19) 

where L = LET represents the linear energy transfer. The graph of Q{L) is illustrated 
in hgure 2. 

Comparisons of the dose and dose equivalent curves are illustrated in the hgures 3 
and 4. The low energy (< 20Mev) neutrons contribute to the dose equivalent next to 
only the primary source protons. HETC-HEDS used the older neutron conversion factors 
to convert from dose to dose equivalent. When the updated conversion factors are used, 
the total dose equivlent from HETC-HEDS increases by approximately a factor of 1.25. 
If this correction is applied to the data in figure 4, then the curve labeled “Corrected 
HETC-HEDS” results and one can see that the three curves are then in good agreement. 

The forward, backward and total neutron fluence curves are illustrated in hgures 5 
through 15. The HZETRN and HETC-HEDS are in excellent agreement since the curves 
are essentially on top of one another over the energy range displayed. The FLUKA results 
are lower over the energy range E < 100 MeV. Also, note that the results from the Monte 
Carlo codes tend to oscillate over the low energy range. This can sometimes be due to 
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lack of convergence of the data, but this may or may not be the actual reason for the 
oscillations. 

In hgures 9-11, HZETRN is differing from both HETC-HEDS and FLUKA over the 
high energy range representation for the backward neutron fluence. This is probably due 
to the HZETRN forward-backward modeling putting too much emphasis on the backward 
component of the fluence. Another possibility for the differences is that the cross sections 
being used for the backward model need modification over this energy range. The rea- 
son for the backward neutron fluence differences will require further investigation of the 
HZETRN code. The HZETRN backward fluence does agree with HETC-HEDS results 
over the low energy range, with the FLUKA results being substantially lower. Note that 
at 30 g/cm^ in the water target, there are no backward neutron fluences. Only outgo- 
ing, forward particles are allowed at this depth. Figures 12 - 15 show agreement of total 
neutron fluence for the three codes over the high energy range E > 100 MeV. Figures 
12 - 15 show differences over the low energy range. The results from FLUKA are lower 
than those of HZETRN and HETC-HEDS. These hgures also show differences between 
HZETRN and FLUKA in the total neutron huence over the mid energy range. 

The proton curves of hgures 16 - 19 are essentially on top of one another, with FLUKA 
illustrating oscillations over the low energy range which can be caused by lack of conver- 
gence or upon the setting of the lower energy limit for the FLUKA code. Note that the 
computed huence spectrum from HZETRN, which at low energies is in a quasi-equilibrium 
state [25], is divided by the stopping power of protons in tissue with a peak displayed 
around 0.1 MeV. This is similar to the classical Bragg peak as described in reference [25]. 
A similar type of statement can be made about the dips illustrated in hgures 20 - 23. 

Figures 20 - 23 are included to illustrate the light ion production produced by the 
HZETRN code. These values can serve as a HZETRN benchmark for any future com- 
parison studies. Also, note that these results, together with the neutron and proton 
production, are major contributors used by HZETRN in the calculations of the dose and 
dose equivalent values. 

Figure 24 is an enlargement of the energy range 1 < E < 100 MeV for the total 
neutron huence from hgure 15. This graph is representative of the existing diherences 
between the results produced by the three computer codes. This region is particularly 
important because this is where the quality factors undergo the most change, as illustrated 
in hgure 2. 

5 Conclusions 

Selected curves have been presented from a comparison study where the transport 
codes HZETRN, HETC-HEDS and FLUKA each analyzed particle transport for a shield 
and target conhguration of 20 g/cm^ Aluminum slab in front of a 30 g/cm^ water slab, 
which is exposed to the February 1956 SPE, as modeled by the Webber spectrum. The 
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graphs presented show that there are many regions where the transport codes are in 
excellent agreement. However, there are some regions where the results from the three 
codes differ. The differences in the dose and dose equivalent curves suggest that a follow 
up study should be performed to determine how ICRP-60 dehnitions of dose and dose 
equivalent are being used by the three codes HZETRN, FLUKA and HETC-HEDS. The 
modeling of the backward neutrons by HZETRN needs to be modihed as the current model 
is placing too much emphasis on backward neutrons at high energies. The cross sections 
used by the three codes HZETRN, FLUKA and HETC-HEDS are different, and how these 
differences affect the computed results requires further study. This research also suggests 
that cross sections used by HZETRN be updated. Future research into a comparison of 
the results from the HZETRN, HETC-HEDS and FLUKA transport codes for the same 
shield/target conhguration, when exposed to the 1977 solar minimum GCR environment, 
will aid in trying to isolate causes for the differences between the codes. Results from 
this study and future comparison studies can provide insight into the HZETRN code 
operation and bring the code to a high standard of performance so that design engineers 
can have conhdence in its abilities to analyze particle transport associated with a variety 
of radiation environments. 
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Figure 1: Webber 1956 integral and differential SPE spectrum with 100 MV rigidity. The 
integral spectrum vertical axis is on the left and the differential spectrum vertical axis is 
on the right. 
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Quality Factor versus LET 



Figure 2: Q{L) as defined by equation (19). 
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Figure 3: Dose versus depth in water after aluminum shield. 



Figure 4: Dose equivalent versus depth in water after aluminum shield. 
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Figure 5: Forward neutron fluence versus energy. Water depth is 0 g/cm^ 
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Figure 6: Forward neutron fluence versus energy. Water depth is 10 g/cm^. 
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Forward neutron fluence versus energy. Water depth is 20 g/cm 
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Figure 8: Forward neutron fluence versus energy. Water depth is 30 g/cm 
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Figure 9: Backward neutron fluence versus energy. Water depth is 0 g/cm^. 
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Figure 10: Backward neutron fluence versus energy. Water depth is 10 g/cm^ 
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Figure 11: Backward neutron fluence versus energy. Water depth is 20 g/cm^. 
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Figure 12: Total neutron fluence versus energy. Water depth is 0 g/cm^. 
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Figure 13: Total neutron fluence versus energy. Water depth is 10 g/cm^. 
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Figure 14: Total neutron fluence versus energy. Water depth is 20 g/cm^. 
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Figure 15: Total neutron fluence versus energy. Water depth is 30 g/cm^. 
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Figure 16: Total proton fluence versus energy. Water depth is 0 g/cm^. 
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Figure 17: Total proton fluence versus energy. Water depth is 10 g/cm^. 
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Figure 18: Total proton fluence versus energy. Water depth is 20 g/cm^. 
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Figure 19: Total proton fluence versus energy. Water depth is 30 g/cm^. 
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Figure 20: HZETRN total fluence versus energy. Water depth is Y (g/cm^). 
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Figure 21: HZETRN total fluence versus energy. Water depth is Y (g/cm^). 
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Figure 22: HZETRN total ^He fluence versus energy. Water depth is Y (g/cm^). 
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Figure 23: HZETRN total ^He fluence versus energy. Water depth is Y (g/cm^). 
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Figure 24: Enlargement of Figure 15. 
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